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Abstract 

Although one can formulate an intuitive notion of instantaneous frequency, generalizing "frequency" as we 
understand it in e.g. the Fourier transform, a rigorous mathematical definition is lacking. In this paper, 
we consider a class of functions composed of waveforms that repeat nearly periodically, and for which the 
instantaneous frequency can be given a rigorous meaning. We show that Synchrosqueezing can be used to 
determine the instantaneous frequency of functions in this class, even if the waveform is not harmonic, thus 
generalizing earlier results for cosine wave functions. We also provide real-life examples and discuss the 
advantages, for these examples, of considering such non-harmonic waveforms. 
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1. Introduction 

The term "instantaneous frequency" is somewhat of an oxymoron. In many cases, researchers who seek to 
decompose signals into different components, and who want to determine the "instantaneous frequency" |15l 
[TB] of each, really seek to solve problems of the following form: 

given a function of the form 

K 

/(i) =^^fe(t)cos(27r0fe(t)), withAk(t),cp'k{t)>{)^t, (1) 

k = l 

compute Ak{t) and (f)'i^{t) and/or describe their properties. 

The Synchrosqueezing transform, a reassignment method [121 [5] originally introduced in the context of 
audio signal analysis [7] , and further analyzed in [51 dni HOI [TS] , provides a way to determine Ak (t) and {t) 

uniquely, up to some pre-assigned accuracy, under some conditions on f{t). 

For some applications, decompositions of the form ([T]) are too restrictive. Consider, for instance, the function 
/ illustrated in the left plot of Figure [T|a); it is a toy example of the type 

fit) = A{t)s{27rm} (2) 

where, as before, A{t) and (f>'{t) vary slowly, and s is now no longer a cosine, but the periodic extension of 
the function in the right plot of Figure[T]ja). It is clear that this f{t) can be decomposed as in ([T]), simply by 
replacing s by its Fourier expansion; however, the representation ([2| is much more efficient because it uses 
fewer terms. Another way in which one could reduce f{t) to an expansion of type with a single term. 
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would be to "absorb" some of the properties of s into a modified phase function m{t) := cos~^ {s{2TT(j>{t))). 
This is the preferred solution in [10], where m{t) is then called the "intra- wave" modulation, indicative of the 
nonlinear nature of the wave process generating f{t). Although this would work for the example in Figure 
[ija), it is not always possible to do this. Consider the function / in the left plot of Figure [TJb), clearly of 
the same type as in Figure [ija); one would expect that the same remarks apply to both examples. However, 
because s has several maxima, it can not be written as cos(m(t)), with m a monotonic mapping on [0, 27r]. A 
single-term expansion of type ([T]) can thus not deal with this example by using only "intra- wave modulation" . 
To nevertheless still reduce f{t) to an expansion of type ([T]) with a single term, one can, instead, attempt to 
"absorb" some of the properties of s into a modified amplitude A. In this particular case, 



sit) 



cos(0.8cos(t)) 



sin(t) 
cos(t) 



sin(0.8cos(i)) - 1.4cos+(t -h 1/6)^ 



cos(i) =: w{t) cos(i). 



(3) 



where cos+(<) — max{0, cos(i)}. So one can write f{t) — cos(27r^(t)), with A{t) = w{2Tr(f>{t))A{t). 
Note that this amplitude A varies much faster than A, masking the only slowly changing wave-pattern of 
f{t). Even this solution is not always applicable. Consider the function / plotted on the left of Figure [ijc). 
Absorbing the "extraneous" extrema of the basic wave-pattern into a special A would lead to an amplitude 
that is no longer always positive, clearly an undesirable trait. 

Although the examples in Figure [T] are just toy examples, similar phenomena can be observed in real-life 
signals, for example, in electrocardiography (ECG) (see Figure [2]). As we shall see below, it is important, 
in these real-life examples, to tease apart the characteristics of the "shape" s{t) from the slow variations 
in and A(t). Rather than insisting on a representation of type ([T]), we are therefore, in this paper, 

interested in decompositions of the type 



K 



(4) 



fe=i 



where we shall give the name "wave-shape function" (or, shorter, "shape function") to the 2tt periodic 
functions Sk, generahzing the cosine functions of ([!]) and 6J, and where we assume that |^fc(i)| and 1 0^(^)1 
are small compared with (/>'j,(t), as in [5]. Functions of type Q can be found in many applications. We 
already mentioned ECG signals; another medical signal of this same form is respiration. 

The additional layer of generality in Q, when compared to ([l]), adds to the complexity of determining 
desirable decompositions of type Q for a (noisy) signal /. Even for decompositions of type ([I]), uniqueness 
is not guaranteed (see e.g. [6]); this absence of unqueness can obviously be only more pronounced for decom- 
positions of type Q . Even when the decomposition is unambiguous, the algorithmic task will necessarily be 
more complex, since not only the A^it) and 0fe(i), but also the Sfe(i) need to be determined, in general. We 
shall concentrate here on the determination of the amplitudes Ak (t) , the instantaneous frequencies (j)'^ (t) and 
the components Ak{t)sk{2TTcj)i~(t)) for the more generalized expansions Q for a restricted class of wave shape 
functions. For more general wave shape functions Sk, the determination of the wave shapes Sk themselves 
will be discussed in a sub-sequential paper. 

The main result of this paper is that, under the same technical conditions as in [B] (the 0'^ have to be suffi- 
ciently separated), we can determine Ak{t) and 0'j,(i) and reconstruct each component via Synchrosqueezing, 
for suitable wave forms Sk- We shall apply this method to ECG and respiration signals, and show the clinical 
potential of the results. In the next section, we first discuss some properties of these signals in more detail. 
Section 3 states and proves the related theorems, which generalize the results provided in [6]. 



2. Two biomedical signals 



The ECG signal can be measured easily and cheaply, and its clinical usefulness is well established. It is 
commonly accepted to model the ECG signal as a current dipole vector undergoing a periodic motion in 



2 




time (second) time (second) 

Figure 1: Toy examples of the type f{t) = A(t)s{2TT4>{t)). The amphtude modulation function and the phase function are 
identical for the three cases: A{t) = %/l + 0.2t and <j>{t) = t + 0.3cos{t). In (a) the shape function s{t) is given by s(t) = 

[cos(0.8cos(t)) - SHiii sin(0.8cos(t))j cos(t); in (b): s{t) = [cos(0.8 cos(t)) - sin(0.8 cos(t)) - 1.4 cos+ (t + 1/6)^] cos(t); 

in (c): s{t) = [cos(1.2 cos(t)) - SiW sin(1.2 cos(t))] cos{i). 

[H]; the recorded ECG signal is then viewed as the orthogonal projection of this dipole vector onto a 
fixed axisj^ The ECG signal, recording the dynamics of the electrical activity of the heart, is a collection of 
periodic oscillating time series, one per channel, each corresponding to an ECG lead [HI IB]- The waveforms 
provide a lot of information about the anatomic or electrophysiological structure of the heart, essential 
in ascertaining certain medical conditions, such as ischemia or atrial fibrillation. On the other hand, the 
variation of the time intervals between sequential heart beats, referred to as "heart rate variability" (HRV) , 
has been shown in the past few decades to be related to more general physiological dynamical processes |13j . 
To understand these different types of dynamics of the physiological system via the recorded ECG signal, it 
thus is beneficial to separate the shape of the oscillation from the variability of the time intervals between 
sequential oscillatory waveforms. 

In this section, the example ECG signals are recorded from a healthy 33- year-old male, with 12 bit resolution 
and a sampling rate of lOOOHz. We show in Figure [2] the lead I and lead II ECG signals, which we denote 
as ECGi(i) and ECGii(t), t £ [0,T]. In the figure, the black curve is the usual ECG signal; the peaks 
indicated by red circles are called the R peaks [8^. One type of variation from one cycle to the next in the 
ECG signal can be clearly tracked by the changing time intervals between consecutive R-peaks. In clinical 
practice, the "heart rate" is given by simply counting the number of beats during a minute, that is, the 
"mean rate" over a given time period. However, there is information hidden inside the HRV, beyond the 
mean rate, that we want to describe quantitatively. To access this, let us first introduce an intuitive definition 
of time-dependent instantaneous heart rate (IHR) , as the inverse of the time interval between the two most 
recent successive heart beats. We refer to this quantity as the intuitive instantaneous heart rate (or intuitive 



^In fact, the situation is a bit more complex. In addition to the periodic motion for every heartbeat, the trajectory in M? of 
the current dipole vector is deformed by other phenomena such as breathing. We shall ignore these effects here. 
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Figure 2: Top: the lead I ECG signal; bottom: the lead II ECG signal, where the red circles mark the locations of the R peaks. 
It is clear that the time interval between consecutive R-peaks is not constant. 

instantaneous frequency) and denote it as IHRi(t), where the subscript i refers to the "intuitive" character 
of this definition. Denoting by tk the location of the R-peak of the A;-th heart beat, k = 1, . . . . IHR^ is 
defined as the piecewise constant function plotted in Figure [3j 



IHR,(t) 



li tk < t < tk+i- 



Note that the definition of IHRi echoes the etymology of "frequency" : counting how frequently a phenomenon 
occurs per unit time. 
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Figure 3: The intuitive instantaneous heart rate. 

To a very good approximation, the ECG signals can be modeled as 

ECG,(i) = A,{t)s,{^E{t)) + W,{t), 

where i stands for the lead (or channel) number (e.g., I or II), and where (jisit) does not depend on £. 
(The subscript E indicates that this is the (t>{t) for the ECG signal; we shall meet other (p below.) The 
function Wi{t) contains all the other components of the signal, including noise and low frequency baseline 
wandering. We show in this paper that the Synchrosqueezing transform can identify cj)'^{t) for such ECG 



^The official code of the Synchrosqueezing transform is available in http://www.math.princeton.edu/-6brevdo/synsq/ 
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signals, with high accuracy, independently of the detailed properties of the wave shape functions si{t), for 
se within a certain class of functions (that contains ECG profiles). Physiologically, the HRV reflects the 
periodic behavior of the current dipole vector; it should thus not depend on how we project this dipole 
vector. In other words, our estimate should be the same, whether we derive it from ECGi(i) or ECGii(t), 
i.e. whether the behavior in time of the wave shape function it modulates is associated with lead I or lead 
II. It turns out that the Synchrosqueezing estimate does indeed exhibit this invariance, as we demonstrate 
numerically here. (A mathematical explanation of this will be given in Section [sj) We denote the estimated 
c/)'^ as IIIRf(t) if it is estimated by computing the Synchrosqueezing transform of ECG^(t). Figures [ij resp. 
[5] plot IHRi(t), resp. IIIRn(i) and show that in both cases IHR^(i) is extremely close to IHR^. In Figure |6] 
IHRi and IHRh are plotted together, illustrating directly that our estimate of (p'^ does not depend on which 
wave shape function was being modulated. (That the two shape functions Si and sn are different is clear 
from Figj2]) Furthermore, this example indicates that our estimate captures the "instantaneous frequency": 
indeed, (f>'^ matches IHR.; with high accuracy. This finding makes it possible to define an "instantaneous 
heart rate" for subjects with much less clearly defined heart beats: even if the IHRi(t) can no longer be 
defined because the tk cannot be defined clearly, we can still estimate IHR£(i). 
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Figure 4: Top: the Synchrosqueezing transform of ECGi(t); bottom: the black curve is IHRj(i) estimated from the Syn- 
chrosqueezing transform of ECGi(t) and the red curve is IHRi(t) shifted down by 0.1. 

Our second example concerns the respiratory signal, i.e., the signal provided by a mechanical recording of the 
breathing process. Physiologically, the morphology of the recorded respiratory signal reflects the anatomical 
structure of our respiratory system [9]. In a first approximation, breathing is a periodic process; a closer look 
reveals that consecutive "breathing peaks" are not quite equally spaced (see Figure [T] below). This breathing 
rate variability (BRV) reflects the physiological dynamics [31 [HI [Ij. An example of a respiratory signal 
Resp(i) is shown in Figure [t] The signal was recorded from a 28-year-old healthy male, with a sampling rate 
of 20 Hz. Like the ECG signal, to a very good approximation, the respiratory signal can be modeled as 

Resp(t) = A{t)s{(l3Rit)) + W{t), 

where the subscript R indicates that this is the phase function (f){t) for the respiratory signal; the function 
W(t) contains again all the other components of the signal, including noise and low frequency baseline 
wandering. 

As in the discussion of the ECG signal, we first define an intuitive notion of time-dependent instantaneous 
breathing rate IRR^ (t) as the inverse of the time interval between the two most recent successive breathing 
cycles. A breathing cycle is defined to be the signal between two consecutive ends of the inspiration; these 
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Figure 5: Top: the Synchrosqueezing transform of ECGn{t); bottom: the blue curve is IHR/j(f) estimated from the Syn- 
chrosqueezing transform of ECGii{t) and the red curve is IHRi(t) shifted down by 0.1. 




Figure 6: The blue curve is IHR/(t) and the black curve is IHRii(t) shifted down by 0.1. 

ends are marked as the red circles in Figure [7] More precisely, we denote the location of the ends of the 
inspiration hy tk, k = 1, M and define IRRi(i) as follows: 

IRRi(i) = \ when tk < t < tk+i- 

The IRRi(t) is plotted as the black piecewise constant curve in FigurejTj Next, we apply the Synchrosqueezing 
transform directly on Resp(i) and get an estimation of its instantaneous frequency, denoted by IRR(t). In 
Figure [7) the IRR(i) is superimposed on Resp(t) and IRRi(t) to demonstrate that the estimated IRR(t) 
captures the notion of instantaneous frequency. Indeed, the spacing of respiration cycles in Resp(i) is 
reflected by IRR(t): closer spacing corresponds to higher IRR(t) values, and wider spacing to lower IRR(t) 
values. In Figure |8] IRR(t) and IRRi(t) are put together for comparison; they are clearly closely related. 

For the sake of convenience, we shall use the acronym SSTIF for the Synchro Squeezing Transform-derived 
Instantaneous Frequency in what follows. Figure [9] shows a different breathing signal, and illustrates a direct 
physiological application of the respiration SSTIF. Sleep is a universally recurring physiological dynamical 
process. It is divided into two broad stages: rapid eye movement (REM) and non-rapid eye movement 
(NREM). Normally, sleep proceeds in cycles, each alternating between REM and NREM, with one cy- 
cle taking about 90 minutes. A clinically acceptable staging of the sleep is determined by reading the 
recorded electroencephalography (EEG) based on the R&K criteria, which were standardized in 1968 by 
Allan Rechtschaffen and Anthony Kales [17] . We take the staging according to this criteria as the gold 
standard, and we demonstrate that the respiration SSTIF can recover the gold standard staging of REM 
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Figure 7: The blue curve is the respiratory signal Rcsp(t); the black piecewise constant curve is the intuitive respiration rate 
IRRi(t); the smooth red curve is the IRR(t) estimated via Synchrosqueezing (or SSTIF). The vertical axis is for IRR(t) and 
IRR,{t). 
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Figure 8: The red curve is the estimated instantaneous frequency IRR(t) of the respiratory signal from the Synchrosqueezing 
transform, and the black piecewise constant curve is the intuitive respiration rate IRRi(i). 

versus non-REM. The result is illustrated in Figure |9] which shows the SSTIF of a respiratory signal recorded 
for about 8 hours with a 16 Hz sampling rate (original signal not shown here), and the corresponding sleep 
stages, denoted by S{t), determined from a simultaneously recorded EEG according to the R&K criteria. 
High correlation is observed when comparing the time intervals with S{t) = 5 with SSTIF. 

The discussion and examples in this Section illustrate the importance of an accurate determination of the 
"instantaneous frequency" of a signal, as captured by our SSTIF notion. In principle, we can determine 
instantaneous frequencies (t)'{t) from representation of type (1) as well as from ([2]) or ([4]). We believe (and 
will argue in the next Section) that, at least for certain signals /(i), modeling / as in ([2]) or ([4|, i.e., in 
terms of "wave shape functions" rather than with cosines, leads to more accurate estimates of the function 
(p'it) via Synchrosqueezing, and probably also by other methods. For such signals, it is thus important to 
separate each component of / modeled as in ([2]) or (|4]) into a "wave shape function" on the one hand, and 
slowly varying amplitude and instantaneous frequency on the other hand. 

Another reason why it is important to separate "shape function" from "instantaneous frequency" is that 
the information hidden in the shape function is important in its own right, and can, for medical signals, 
be used for clinically quite different diagnoses. In other words, this separation allows us to tease apart two 
types of information that are commonly mixed-up. For example, reading EGG signals to diagnose cardiac 
disease in clinical practice amounts to evaluating the morphology of each heart beat, and this is given by the 



^Further study and finding of the sleep cycle is beyond our scope here; a more detailed study will be presented in a later 
paper. 
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Figure 9: The blue curve is the lOxSSTIF+4, and the red piecewise function is the sleep stage determined from the simulta- 
neously recorded EEG by the R&K criteria. The REM stage period is emphasized in both signals. 



wave shape function. A well known example is the typical "ST elevation" in myocardial infarction patients. 
Another example is the qualitatively different spectral behavior of the morphology of the ECG shape function 
indicative of myocardial ischemia: the spectral analysis of ECG wave shape functions from dogs revealed 
a shift from high- to low-frequency ranges in ischemia cases [14] . Similar phenomena have been associated 
with balloon inflations during percutaneous transcatheter angioplasty in CAD patients We expect that 
a cleaner identification of the SSTIF will also lead to a separation of the shape function particularities with 
greater sensitivity, which in turn will be useful for clinical diagnoses linked to those particularities. A detailed 
discussion of this will have to wait for a subsequent paper; here we concentrate on the identification of the 
instantaneous frequency. 



3. Theorem 



We start by introducing some notations and conventions. Suppose / is a tempered distribution. The Fourier 

transforms of / is denoted by / and satisfies the normalization e~'^^^ — e~^^ . In the following, we denote 
Wf{a,b) to be the continuous wavelet transform of / with the mother wavelet ip which is in the Schwartz 
space [5]. For simplicity, we shall assume supp ip C [1 — A, 1 -I- A], where < A <C 1; in practice, our results 
seem to hold under much less stringent conditions. 

As is customary, we say that r is a period for the function (/? : M — >■ C, or that ip is r— periodic, if, for all 
<: G E, and all fc S Z, (p{t + kr) — cp{t). We shall designate T to be the Period for (p, if 

Vi e M , Vfc e Z : ip{t + kT)= (p{t) and T = inf ({ r ; r > and r is a period for ip }) . 

Next, we define a special class of wave shape functions that are dominated by one particular Fourier mode. 
More precisely. 

Definition 3.1 (Analytic shape function class S^'^'^). Fix 5 > 0, D S N and 6 > Q. The analytic shape 
function class S^'^'^ is the subset of C^'"(T), where T is the 1-dim torus and a > 1/2 (the 27r— periodic 
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a-Holder continuously differentiable functions), consisting of analytic functions s with mean 0, i.e. s{k) = 
for all fc < 0, and unit L^— norm for which all the Fourier modes ?(fc), k ^ I are dominated by the product 
of S and the first mode coefficient, i.e. 



and 



Vfc e N, with k^l, \s{k)\ < S |s(l)| 



^ \ns{n)\ < 9 

n>D 



(5) 



(6) 



Remark. (1) The parameters 6, D and 9 are used to characterize the "shape" of the shape function. Condition 
([5]) says that the "base" frequency s(l) can not be zero; we use this condition to estimate the "instantaneous 
frequency" as we should see below. Condition ^ says that "essentially" the shape does not oscillate too 
fast. (2) In practice, we shall consider real-valued shape function s, which are not analytic. However, we 
can use the standard trick of viewing them as the real part of an analytic function, i.e., s = 5fts, where 
s{k) = s{k) if A; > 0, s{k) = if fc < 0. Because the Synchrosqueezing transform (to be defined soon) uses 
the positive frequency domain only, s and s can be used interchangeably in our analysis. Thus, we can define 
a real shape function class S^'^'^ of parameter S, D, and 9 to contain the functions s : T — ^ K so that s = Sfis 
for some s £ S^'^'^ . Since the analysis of the functions with the real shape functions is the same as that 
of the analytic shape functions, in the following we focus on analytic shape functions in our analysis, for 
notational convenience, but when we demonstrate numerical results, we use real shape functions, which is 
easier for visualization. 



For example, e** is a shape function, which is widely used in Fourier analysis. Indeed, e**(n) = 1 when ri = 1 
and when n ^ 1 so that e** S 5°'^'°. We demonstrate some examples in Figure 10 the shape function 
su{t) for the ECG lead II signal satisfies Su{t) e 5*"°'^ where 6 « 3.4, L» « 40 and 6^ < 1. (Note that the 
shape function Su{t) depends on the subject; different ECG signals may belong to different S^'^'^; the value 
i5 = 3.4 is one of the larger we have encountered - more often S can be picked less than 1.) 




1 S 1 A i S 



Figure 10: Three different shape functions si{t), S2{t) and sii{t) (upper) and their spectrum (below). Left: si{t) is equal to 
normalizing |^cos(0.5 cos(t)) — sin(0.5 cos(t))j cos(t) to be of unit norm; middle: S2(t) is equal to normalizing Jsl to be 

of unit norm; right: su is a shape function for a particular lead II ECG signal. 



Next, consider the following class of functions. 
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Definition 3.2 (Intrinsic Mode Functions class C^'^). For fixed choices of e, (5, D > 0, e ^ 1, the space 
Cf'^ of Intrinsic Mode Functions (IMFs) consists of functions / : M -> C, / G Ci(M) n L°°(M) having the 
form 

f{t) = A{t)s{2irc^[t)), (7) 
where s € 5*'^'^, such that A and (j) satisfy the following conditions: 

A G Ci(M) n L°°(M), 4) e C^{R), 

inf^(t)>0, inf0'(t)>O, sup (/)'(t) < oo 

*GR teR igH 

< e\(l)'{t)\, \(f>"{t)\ < e|0'(t)| for aU t G M, 

M" := sup 10" (01 < oo. 
teK 



We then consider the function class C^'^, defined as follows. 

Definition 3.3 (Superpositions of IMFs). The space Cf'^ of superpositions of IMFs consists of functions / 
having the form 

k=l 

for some finite K > and fk{t) — Ak{t)sk{2n(j)k{t)) G C^'^ such that (f>k satisfy 

Mt) > <Pk-iit) and </)Ut) " <^'fc_i (0 > d[<P',it) + 0Li W]- 



Note that Cf'^ and S^'^'^ are not vector spaces and are not the most general possible definition. We use 
these definitions for the purpose of simplifying the proof and clarifying the main idea of this work. 

In the following we prove that the Synchrosqueezing transform discussed in ^ allows us to estimate the 
instantaneous frequency of each component of the functions in C^'^ with high accuracy when all the wave 
shape functions are close to the exponential function, i.e. when 6 is sufficiently small. Moreover, under some 
constraints on the IF of each component, we are able to reconstruct each component. Most of the estimates 
needed to prove these results follow from Section 3 in 6]; we shall spell out the details only when care has 
to be taken with extra terms involving the possibly more complex shape of the wave shape functions. 

Definition 3.4 (Instantaneous frequency information function). Let / G L°°(M). The instantaneous fre- 
quency information function of / is defined by 

4"^^ i«',«..i.)i>o 

loo \Wf{a,b)\=0 



Note that this definition makes sense since Wf{a, b) G C°°(M+ x M). 

Remark. In practice, the determination of those (a, 6)-pairs for which Wf{a, 6) = is rather unstable when 
/ has been contaminated by noise. For this reason, it is useful to consider a threshold for |W/(a, b)\, below 
which u)f{a, b) is not defined. The purpose of this function is to record the information of the instantaneous 
frequency, based on which the reassignment will be performed. 

Then, we consider the following definition of the Synchrosqueezing based on the wavelet transform: 
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Definition 3.5 (Synchrosqueezing transform). For / G L°°(M), the Synchrosqueezing transform with reso- 
lution a > and threshold 7 > is defined by 



JA-,j{b) a \ a J 



(8) 



where 6 e M, ^ e M+, A^j{b) := {a G M+; \Wf{a,b)\ > 7} and h{t) = -^e~ 



Note that ujf{a,b) and S'"'''(6, ^) can be defined for any L°° function /. However, when / e C^'d 1 the next 
Theorem tells us that the Synchrosqueezing transform provides an accurate estimation of the instantaneous 
frequency and allows the reconstruction of each component. 

Theorem 3.6. Let f{t) = J2k=i Au{t)sk{2'K(t>k{t)) G Cf'^. Suppose A < d/{l + d) and 



l-d 1 



n<P'^{ty n<j,'^{t) 



(9) 



for aU 1 < fc < A' and t e M. Set := / i(j{z)z-'^dz, Zk,n := {(a, 6) : \a.n(|)'^.{b) - 1| < A} and e := e^/^. 
Then, provided e is sufficiently small, the following hold: 



|W^/(a, b)\ > € only when (a, b) G Z/c,n for some fc € {1, . . . , and n e {1, 
For each pair (a, b) G Zfc^„ for which |T4/'/(a, b)\ >e, we have 

\ujf{a,b)~n^'^{b)\<Z 

Moreover, for each k S {1, . . . , K} and all 6 G M, 



,D}. 



lim 7^Tl 



a-j-O 



where C = Afc(6)e^ + 4 



1/2 /^/,M\l/2' 



< Ce 



Notice that the special case when if = 1 is important since the EGG signal and the respiration signal 
are both signals of this kind. The results of estimating the instantaneous frequency via Synchrosqueezing 
transform have been shown in Section 2. 

When the wave shape functions of all the components of f{t) are close to the imaginary exponential function, 
the following corollary immediately follows by taking D = 1 in Theorem |3.6[ 



Corollary 3.6.1. Let f{t) = T,k=i Ak{t)sk{2T:(j3k{t)) £ C"'^. Suppose A < rf/(l + and (5 < e. We use the 
same notations 7?,^, Z^^n and e^as in Theorem 3.6 Then, provided e is sufficiently small, the following hold: 

• I W/(a, 6)1 > ? only when (a, b) G Z^^i for some fc G {1, . . . , K}. 

• For each pair (a, b) G Zk^i for which |M^/(a, b)\ > e, we have 

|a./(a,6)-0U6)| <?. 

• For each k e {1, . . . ,K} and all 6 G M, 



(5,1 



lim 7^Tl 



where C = Ak{b)'P + 4 



{eiC-0'fc(b)l<?} 

1/2 /w//M\l/2' 



Sf'{b,Odq - Ak(b)sk{2TTMb)) 



11 



Notice that in Corollary 3.6.1 the wave shape functions of all components are close to e** while the separation 
of the "base frequencies", that is, (j>'i^{t), are not further limited; on the other hand, the wave shape functions 
of each component in Theorem 3.6 is more general than the cosine function but the behaviors of (j)'i^(t) need 
to be further limited by Setting Sk{t) = e** for all fc = 1, . . . , if, we fully recover Theorem 3.3 from 

-,S,D 



Condition (joj) can be shghtly relaxed as follows. Suppose f{t) — J2k=i ^kit)sk{'2n(j)kit)) e 
Sk G S^'^'"'^ and Dk < D maybe different. Then the following condition suffices for the proof: 



1+d 



D, 



l-d 1 



where 



(10) 



Since the proof is entirely analogous to what we show below, involving longer computations without making 
a conceptual difference, we focus here on the condition ^ only. 

The proof of the Theorem is similar to the one carried out in [6, Theorem 3.3], except for some estimates 
related to the wave shape function, which are given in detail below. In the statement and proof of all the 



Lemmas, we shall always assume that the conditions of Theorem 3.6 are satisfied without repeating them, 
unless stated otherwise. 

First of all, we have the following estimates bounding the growth of Af^ [t) and (j)'^. (t) in terms of the value of 
(p'kit)- The proof is the same as that of Estimate 3.4 in so we omit it. 

Lemma 3.7. For each k E {1, . . . , K}, we have 



\A,{t + u) ~ Ak{t)\ <e\u\( \cl,',{t)\ + -A^'H 



and \Mt + ^) - Mt)\ < \Mt)\ + 7;M'^\u\ 



The next lemma concerns the "dyadic separation" condition in the definition of C^'^ . It implies that for 
every (a, 6)-pair, at most one component of the signal "comes into play" . 

Lemma 3.8. For any pair (a, 5) such that a > jj^p-^, there can be at most one / S {1,...,K} and 
ne{l,...,D} for which |an0J(fe) - 1| < A. 



Proof. Suppose ^ holds. We rewrite f{t) as 

K K D K oo 

fc=l fc=ln=l fc=lri=D + l 

where the pointwise convergence holds since G C^^". By condition we know for all t G 



(11) 



l-d l+d 



l-d l+d 



mcfy'Aty m(t)'At) 



(12) 



when / 7^ J or n 7^ m, where n, m < D. Next we permute D x K components {n(j)'i^{b)}n=i^...^D,k=i,....K in 
the ascending way and index (j)'i{b) by 1 and D(t)'j^{b) hy D x K. Thus (11) can be further rewritten as 



DxK 



K oo 



^ ii(i)e*2'''^'(*) +E E ^fe(<)sfc(n)e'2™'^''(*\ 



1=1 



k=l n=D+l 
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where Ai{t) = si{l)Ai{t), . . . , AoxKit) = SK{D)AK{t) and ViW = <hjt), ■ ■ ■ ,i'DxK{t) = D(t,K{t). Take 

l-A _ 1-A 



j,l e {l,...,D X K} SO that / > j and fix a > ^^r^ = wt3^- l^^l* '^^'^ 

^P'i{h)-xP',{b)>d[^r(b)+^:[{b)] 



(13) 



since < ^|r^- Now suppose 107/1^(6) — 1| < A and \atl)[{b) — 1| < A hold simultaneously. A direct 

calculation leads to 

/Tr /J 



and 



V';(6) + V''(6)> 



which with ( 13 ) gives 



(1-A) + (1-A) 2(1 -A) 



A > d(l - A) 



contradicting the condition A < 



+d- 



□ 



The following Lemma describes the result of applying the continuous wavelet transform to / G C^^'d- 
similar to Estimate 3.5 in [B] except for some extra terms due to the more general form of the wave shape 
functions; when Sk{t) = e** for all fc G {1, . . . , K}, we recover the statement and proof of Estimate 3.5 in [B]. 

Lemma 3.9. For Z e {1, . . . , A'}, ri e {1, . . . , £)} and (a, h) E Zi n, we have 



Wf{a, b) - A,(6)si(n)e"^"'^'('')v^7/. {an(l)[{b)) < ea^'^Ki{a, b) 



(14) 



where 



Ai(a,&) = ^{||sfc||^ U;.(6)a/i + 



fc=i 



2 



-7r^fc(&) 



3 



, }+/oEA,.(6) 
^ fe=i 



and li ~ |xp 1-0(2;) |dx. 

Remark. We observe from this Lemma that when ^ holds, the larger the I is, the larger the Ai(a, 6) in 
Lemma 



3.9 



Indeed, take / — J2k=i fk G ^t'd ^^'^ -f^ ^'^^ some I € {1, . . . , K} for example. The larger 
I is, the more (f>'i.{b) will be smaller than 4>'i{b) so that more il){a(t>'^.{b)n), n > D become (possibly) positive. 
Consider a such that a(f>'i{b) e [1 - A, 1 + A]. For fee {1, . . .,K}, k < I, i){a(j)'^[b)n) ^ holds only if 



n.<^',(b)(l-A) 
<t>[{b) 



<l + Aorl-A< 



«0;,(b)(i+A) 



, or equivalently 



i-A0;(b) 

1 + Af,(6) 



< n < 



1-A0;(6)- 



(15) 



Clearly when the I is large, '0(a(/)'j,(6)n) 7^ for more combinations of k and n. On the other hand, if 
Skit) 7^ e** for some k < I, due to the existence of the nonzero high Fourier modes of the shape function 
Skit), that is, Sk{n) ^ for n > 1, the nonzero ^{a<f>'^.{b)n) will survive for those n for which Sk{n) ^ 0. 
Thus, the continuous wavelet coefficients around n = 1, . . . will be "contaminated". Putting these 
two effects together, we expect that the larger / is and the more non-zero Sfc(n), n > D and k < I are, the 
larger the Ai(a, 6) will be. Thus, the higher the n > 1 is, the more blurring is around the band Zk.n on 
the time-frequency plane, which depend on the wave shape functions. It is for these reasons that we shall 
evaluate the IF from the lowest mode Sfc(l) of the Sk- We shall illustrate this below, in examples. 
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Proof. Since / G Cf'^f C and the mother wavelet ■(/'(t) is a Schwartz function, Wf{a, b) is well defined. 

We work out the estimation step by step. First, we evaluate the following integration for fc e {1, . . . , K} and 

(a, 6) e M+ X M: 



dt 



(16) 



-dx 



-dit 



^i2Trn{(t>k{b) — u) 



-du 



where the first two equalities comes from the change of variable and the third equality comes from the fact 
that sk e Ci'"(T). 



Next, by applying Lemma |3.7[ we can evaluate the difference between the continuous wavelet transform of 
Afe(6)sfe(27r0fc(i)) and Ak{b)sk(2^{d^k{b) - b<j3'k{b)) + 2^0'fc(6)t): 



At{b)sk{2TT(j)k{t))^i} ^ ^ 



/a \ a 



dt 



A 



k{b)sk(27T{Mb) - bMb)) + 27r^'fe(&)i) 



1 , /t-b 



a \ a 



dt 



Akib)sk 



27rlMb)^b<P',{b)+<P',{b)t+ [<P',{b + u)-<t>Ub)]du 



t-b 



\/a \ a 



dt 



Ak{b)sk 2TT{Mb)-bMb)+Mb)t) 



1 , /t-b 



/a \ a 



dt 



= Akib) 

< Ak{b) 

< 2Tr Ak{b) 

< 2TTeAk{b) 

< 2iT€Ak{b) 



1 , (t-b 

--i' 



'a \ a 



dt 



Lne 



51 

nez 



'27™ /o ''[0fc(ft+«)-</''fc(h)]du _ 



t-6 



a 



dt 



h^t-b\^\c^',(b)\ + \\t-b\Hl'^ 



la'/^I,\^',ib)\ + \M'la^'^h 



1 



dt 



t-b 



dt 



where the first equality comes from applying Taylor's expansion to ipkit) and the fifth inequality conies from 
Note that X^nez 1^1 l^fcl")! < °° since Sfc G C^'" with a > 1/2. 



Lemma 



3.7 
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Third, we approximate Ak{t) by Lemma 3.7 

K 



Wf{a,b)~y] [ Akib)ski2nMt))^^J ( — ]dt 
^.^i-Js. y/a \ a J 

K 

< / \Ak{b) - A,{t)\\s,i2nMm^ 



(17) 



1 — 1 ''l^ V 



K 



t~b 



dt 



< e 



where the second inequality comes from Lemma |3.7| 



When l£ {1,..., K}, n £ {1, . . . , D} and (a, 6) e M+ x M such that \an(j)[{b) - 1| < A, ([16|) and ^ together 
lead to the result: 



< 



\Wf {a, b) - Ai {b)si (n)e*2.™0, (b) | 

K K 

J2 Afc(5)4(n)e^2^"^'=WV^V;(an</.;(6)) + ^ ^ Ak{b)sk{j)e'^'''^-^''^ V^^ajcj^'kib)) 

k=l,k^l k=l jetij^n 

^ / 1 



(18) 



fc=i 



K 



-2neY, Mb) 



fe=i 



ia^/2/2|0l.(6)| + ^A/^V2j3 



By the definition of IMTs, the second term and the third term in ( 18 1 are bounded by 

f K 



,fc=l 



Moreover, when ([9]) holds, due to Lemma 3.8 and the assumption of the first term in (181 is bounded by 



K 



K 



Y Afe(6)s,(n)e^2.n0.(6)^^(^^^,^(^))^^ ^ Ak{b)SkU)e'^''"*"'^''^V^rP{aj^'k{b)) 

k=l,k^l k=l jefi,j^n 



K 



- J2 J2 Ak{b)skU)e''''"f"'^''^VarPiamb)) 

k=l jefij^n 

l-l l-l 

k=l jeN,j>D k=l 



(19) 



where the first equahty comes from the fact that (a, 6) G Z;_„ and the assumption ([9]), the first inequafity 
holds since \^p{aj(j)'f.{b))\ < Iq and Lemma 3.8 and the second inequahty holds since / G C^'^ ■ To be more 



precise in the first inequahty, when k = I, ip(aj4)'i{b)) ^ holds only for j = n; when fc > /, no € N leads 
to ip{aj(j)'i.{b)) ^ 0; when k < I, ^{aj<j)',,{b)) 7^ is possible only when j > D. The proof is thus done. □ 
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In the next Lemma, we show that by differentiating the continuous wavelet transform, the instantaneous 
frequency and its multiples pop out. Though the instantaneous frequency is mixed up with other 

quantities, we will show that, when handled properly, this gives us an estimate of the instantaneous frequency 
with high accuracy. This Lemma is analogous to Lemma 3.9 in 



Lemma 3.10. For ^ e {1, . . . , K}, n € {1, . . . , D} and (a, b) £ Zi n, we have 



'idbWfia, b) - 2^TAl{b)sl{n)e'^''"''''^''^n(f>[{b)^/^i: {an(j)[{b)) < ea^/^A2{a, b), 



(20) 



where 
A2(a,5) 



Y.\^\\U3)\ 



Proof. The proof follows the same lines as that for Lemma |3]9j Under the same conditions, we can evaluate 
the following approximations. First, 



n6N 



t-b 



dt 



(21) 



Second, 



(22) 



A,ib)ski2niMb) - bq^',{b)) + 27r</>',(6)i)^V'' 



dt 



< 2TTAk{b) f 

JR 

< 2Tre Ak{b) 



nez 



b'^{b + u)~Mb)]du 



1 



t-b 



dt 



where J^nez, l*^! I'^fe("')l since Sk € C^'" and a > 1/2. And third. 



< e 



-df,Wfia,b)-y] f Ak{b)sk{2nMt))472'^' 



t-b 
a 



dt 



(23) 



k=l 



When I € {1,. . . ,K}, n e {1,. . . ,D} and (a, b) € 
lead to the result: 

i2TTn<pi(b) ,/ ( 



such that \an(t)[{b) - 1| < A, (21)-(23l together 



idbWf{a,b) - 27rAj(6)s,(n)e*27r„0K6)„0j(5)y^^(a„0j(fo))| 
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K 



K 



k=l,k^l k=l jefij^n 



K 



k=l 



Y.\^\\h{j)\ 



By the same argument as that for Lemma 3.9 we get (20 1. 



□ 



The same remark for Lemma |3.9| holds for Lemma |3.10[ That is, the larger the difference between the 
wave shape functions and the cosine function, the larger the possible negative effect on the precision of the 
estimate. The following Lemma clarifies the role of the function u!f{a,b). Indeed, it states that ujf{a,b) 
provides the information of the instantaneous frequency of each component. Since the proof of this Lemma 
is the same as that of Estimate 3.8 in [B], we shall skip it. 

Lemma 3.11. For I e {1, . . . , K}, n G {1, . . . , D} and (a, b) E Zi ^ so that |W/(a, 6)| >e, we have 

|w/(a, b) - n<j)'i{b)\ < e^/^a^/^{Ain(j)'i{b) + A2/2^). 



It is now clear that with appropriate restrictions on e, the first two claims in Theorem |3 . 6| are proved. Indeed, 
note that by the definition of Zk,n and the uniform lower and upper bounds on we know that the 

Zk^n are uniformly bounded. Then, supposing that |9| holds, due 



values of a for which (a, b) € U^^^ ^k=i 



to the uniform boundedness of Ak{b) and (j)'i^(b), it follows that Ai(a, b) and A2(a, b) are uniformly bounded 
in U^^i Zk^n as well. Thus, for (a, b) £ ^n=i ^k=i ^k.n, there exists e > so that 



(24) 



which leads to ea^/^Ai (a, b) < e. If fvuther we impose the condition that for alH = 1, . . . , X and n — 1, . . . , D, 

e < a-3/2(Ai„0;(6) + A2/27r)-3, (25) 
then e^/'^a^/^(Ain05(&) + A2/27r) < ?. Thus the first two claims in Theorem 



3.6 



hold. 



The final Lemma concerns the reconstruction of each component; when Sk{t) — e** for all k — 1, ...,-ftr, we 
recover Lemma 3.9 in [B]. 

la 3.12. Suppose 
, K} under consideration. 



Lemma 3.12. Suppose that both (24) and (25) are satisfied, and that, in addition, for all b and k £ 
{1, 



e < min 



d^q^[{b) + q>'^{b)r 27d^cj,[{br 



(26) 



Then for any b G 



where C = Ak{b)'P + 4 



1-A 



U°^i{5:|5-«0^(6)|<?} 
1/2 /^',M\l/2' 



5?'^(6,C)dn - Akib)ski2TrMb)) 



4>'M 

l+A 



<Ak{b) + A 



1-A 



1/2 



0'fc(b) 
l+A 



< c?. 



1/2- 
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Proof. Fix & g M. As a function a a, Wf{a,b) e C°° {A-^ f{b)), so by the definition of the wavelet Syn- 
chrosqueezing transform, as a function of S'?''(5, ^) G C°°(M). Thus, we have 



lim / S'?''(6,e)d^ 
hm / / 

a-^/^Wf{a,b) 



(27) 



Wf{a,b)-h 
a 



\^-ujf{a,b)\ \ 



hm 



hm a-^/^Wf{a,b) 



l^ n^~u;f{a,b)\ 



a-i-O 



1 



u?=i{C:|?-"0',(6)|<?} a 

Wf{a,b)a-^/'^da 



\^-cvf{a,b)\ 



Made 
d^da 
d^da 



A^j{b)n[u^^i{a:\Ljf{a.b)~n^'^{b)\<7}) 

where we have used Fubini's theorem for the second equahty, the Dominant Convergence theorem for the 
third equahty, and the approximation of identity for the fourth equahty. Indeed, the integrand on the 
third hne is bounded by a~^^'^\Wf{a,b)\ e L^{A-^j{b)) and converges almost everywhere to a~^^^Wf{a,b) if 
|w/(a, b) — n(j)'f^(b)\ < e^for some n £ {1, . . . , D}, and to zero otherwise. 

We now claim that there is only one I e {1, . . . , K} and n G {1, . . . , D} for which |an(/)J(6) — 1| < A and 
\ujf{a, b) — n(j)'i{b)\ < ?hold simultaneously. In fact, if there exist I' E {1, . . . , K} and n' e {1, . . . , D} so that 
\an' (j>'i, {b) — 1| < A holds, where I ^ V or n ^ n', we get 

|c./(a,6) > - n>;,(6)| - \u:f{a,b)-n<ly\{b)\ > d[ncf,[{b) + n'<l>[,{b)] ~I, 



where the second inequality comes from ([13j) and Lemma 3.11 Notice that 

d[ncj,[{b)+n'cj,[,{b)] > dmin{[0'i(6) + 0^(6)],3(/.'i(6)} > 2?, 



(28) 



where the first inequality holds since n(l)[{b) > 4>'i{b), n'ip'i,{b) > <j>'i{b) and <j>2{b) might be larger than 2(t)\{b) 
and the second inequality comes from ( 26 ) . Thus we conclude that 

\LOj{a,b)-n'<P'i,{b)\>I, 



which is absurd. Next, from Lemma 3.10 and (24) we know that |W/(a, 6)| > e only when |a7i0J(6) — 1| < A 
for some ^ S {1, . . . , K} and n e {1, . . . , Z)}. Hence we know 

A~,j{b) n (U^^i{a : |w/(a,5) - n^'^{b)\ <Z\) = Arj{b) n (U^^i{a : \ancj>'^{b) - 1| < A}) 



and the right hand side of ( 27 1 becomes 



Wf{a,b)a-'^l'^Aa 



A-,j(b)n{u^^,{a:\an^'^{b)-l\<A}) 

Wf{a, b)a-^l'^Aa- 

U°=i{a:|a«0;^(;))-1|<A} 



(u°=i{|a«0;^(b)-l|<A})\Aj.^(b) 



Wf{a,b)a^^''^da 
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Thus we obtain 



< 



n. 



< 



uS^^{a:\an4'^(b)~l\<A} 

{uS^,{\an^i{b)-l\<A})\A,j{b) 
D 



Wf{a, 6)a-3/2da^ - Afe(6)sfe(27r(/.fc(6)) 
Wf{a,b)a-^/^da 

„=1 \J|an0;^(b)-l|<A / 



J|a„</.'J6)-1|<A J(u°^i 



'(u°^a|an0'Jfc)-l|<A})\Aj,j(b) 



Wf{a,b)a-^/^da 



< 



D 



E / -"'''da 

„=l"'k»0'fcW-l|<A 



(29) 



where the second inequahty comes from Lemma |3.9| and |24| and the third inequahty comes from Lemma |3.9| 
and the fact that 

\an<P'^{b)~l\<A J|C-1|<A 



The first term in ( 29 1 is bounded by 

CXD 

Ak{b) 4(n)e^^-"^'=(^) 



n=D+l 



n=D + l 



D 



The second term in (29) can be worked out exphcitly: 

„=i-'l™*L('')-i| 

We thus conclude 



1|<A 



1- A 



1 + A 



hm7^-l / Sj-'{b,0d^-Ak{b)sk{2TTMb)) 



< eAk{b)+Ae 



1- A 



1/2 



1 + A 



1/2- 



□ 



Combining Lemmas |3.9|3.12| completes the proof of Theorem |3.6| 

Remark. This Theorem suggests that if we know a priori that the components of the signal have wave shape 
function close to an imaginary exponential (or if its real part has components close to a cosine function), 
then focusing on the band around the instantaneous frequency (Ij'kib) and its multiples allows us to recover 
the signal. 

Now we demonstrate some numerical results of applying the Synchrosqueezing transform to analyze the 
functions in C^'^ . Take the phase functions <j)i{t) = 1.5t + 0.2cos(t + 1) and (f>2{t) — 4.5(i + 0.2cos(i)) 
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and the amplitude modulation functions Ai{t) = 1 + 0.1sin(t^-^) and A2{t) — ^1 + cos(t). Consider the 
following two functions for comparison: 



and 



/2(t) 



Ml 

3.5 



cos(27r,^i(i)) + ^2(i)si(27r02(t)), 



where si{t) and sii(t) are the shape functions demonstrated in Figure [To| Clearly, /i and /2 are both 
composed of one component with instantaneous frequency 1.5 — 0.2sin(t + 1) and one component with in- 
stantaneous frequency 4.5—0.9 sin(t) but with different wave shape functions in the low frequency component. 



As is shown in Figure 10 the coefficients of the high Fourier modes of sn are quite significant, whereas those 
cos(i) are all 0. In this case, since 02 is roughly 3 times it is the coefficient of the third Fourier 



mode of su that matters. Thus, according to Theorem 3.6 we expect to have worse (f>2 estimation from fi, 
which is borne out by Figure [TT] 
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Figure 11: Left top; the /i signal; left middle: the Synchrosqueezing transform of /i, where the y-axis is demonstrated in the 
log scale; left bottom: the blue curve is 4''2{^)y while the red curve is the estimated (pi^', right top: the /2 signal; right middle: 
the Synchrosqueezing transform of /2, where the y-axis is demonstrated in the log scale; right bottom: the blue curve is </>2(t), 
while the red curve is the estimated 02 . It is clear that the estimation of ipi^ is worse in /i. Notice that the dominant curve 
around (p'^ in the left middle figure is not very marked. This is caused by the small coefficient of the first Fourier mode of sn. 



Next we show the component reconstruction results stated in Theorem 3.6 Take the phase functions (jjsit) = 
t>4{t) = 5t + 0.2cos(<)^ and the wave shape functions S3{t) and S4{t). The Fourier 



4t + 0.1cos(t + 1) and 
modes of S3(i) and S4{t) are illustrated in Figure 12 Consider the function 



where 



and 



/3W = /3,lW + /3,2W, 



fsMt) = A2{t)s4{2TTMt))- 

Since /3,i(0 G and hMt) G know f^it) £ C^^j'^ where e 

take A = 0.1 in the simulation, and it is clear that the intervals 



l-A l+A 



2/25 and d « 1/9. We 
for n — 1 , . . . , 4 and 



i-A 



l+A 



for m 



1, 



3 do not overlap, so the condition (|9|) is satisfied, and hence we expect to 



be able to reconstruct fj.iit) ^-nd /3,2(i) from via Synchrosqueezing transform. The extraction result is 



20 



shown in the left panel of Figure 14 We also demonstrate the robustness of Synchrosqueezing transform in 
this case. Consider 

/4 = /3 + 'tW, 



where W is the white noise with variance 1 and a = ^ var(/3)10i/'^, that is, we add — 2dB white noise to /a 
if the signal-to-noise ratio is defined by 



signal to noise ratio (dB) = 10 log^ 



var(/3) 



We first demonstrate the Synchrosqueezing transform of /j and the reassigned Morlet scalogram |2] in 



Figure 13 Then we extract the components /a^i and /3_2 by applying Synchrosqueezing transform on f^. 



The result is shown in the right panel of Figure [T4| Since reconstructing the components from the reassigned 
Morlet scalogram is not guaranteed, we only demonstrate the reconstruction results of the Synchrosqueezing 
transform. 






Figure 12: Top row: the S3 (left) and S4 (right) wave shape function; bottom row: the first 20 Fourier modes of S3 (left) and 
S4 (right). Note that 53 (n) and Si(n) are not zero but very small when n > 4. 



21 




Figure 13: Left: the Synchrosqueezing transform of f^; right: the reassigned Morlet scalogram. Notice that the performance of 
the Synchrosqueezing transform is equivalent to the reassigned Morlet scalogram for the purpose of extracting the instantaneous 
frequency. 




3456789 19 3456789 10 



Figure 14: Left top: the /a signal; left middle: the blue curve is Ai{t)s:j{2TTip3{t)) , while the red curve is the reconstructed 
^l(i)^3(27r</>3(t)), which is shifted down by 0.3 to improve the readability; left bottom: the blue curve is A2{t)s4{2iT(p4(t)), 
while the red curve is the reconstructed A2{t)s4{2iT<j>4,{t)) , which is shifted down by 0.3 to improve the readability; right top: 
the /4 signal; right middle: the blue curve is j4i(t)s3(27r03(i)), while the red curve is the reconstructed Ai{t)s3{2TT<f>3{t)), which 
is shifted down by 0.3 to improve the readability; right bottom: the blue curve is j42(t)s4(27r<?!>4(t)), while the red curve is the 
reconstructed A2{t)s4{2n(f>4{t)), which is shifted down by 0.3 to improve the readability. 
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